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Abstract 

A system of partial differential equations describing the spatial 
oscillations of an Euler-Bernoulli beam with a tip mass is considered. 
The linear system considered is actuated by two independent controls 
and separated into a pair of differential equations in a Hilbert space. 
A feedback control ensuring strong stability of the equilibrium in the 
sense of Lyapunov is proposed. The proof of the main result is based 
on the theory of strongly continuous semigroups. 

1 Introduction 

Dynamical models of flexible-link robot manipulators are generally described 
by a set of coupled ordinary and partial differential equations, that gives rise 
to series of mathematical control problems in inflnite dimensional spaces [31 
in El El [12] ■ However, flnite dimensional approximate models obtained by the 
assumed modes and flnite elements methods are used more frequently for 
solving the motion planning and stabilization problems [Il[lT]. It should be 
emphasized that the majority of publications in this area is concentrated on 
planar manipulator models with a free end. To study spatial manipulators 
with a tip mass, the mathematical model that describes the motion of a 
multi-link manipulator under the action of gravity and controls (torques and 
forces) was proposed in [T3] . 

The goal of this paper is to study the stabilization problem of the control 
system derived in [13] for the particular case of a manipulator with one 
flexible link. 
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2 Equations of Motion 



A mechanical system consisting of n Euler- Bernoulli beams and a rigid body 
as a load was introduced in [13]. In this paper, we assume that n = 1 and 
neglect controlled rotations of the load {ifj = and c = in the notations 
of [13]). Thus, the beam deflection at time t is defined by functions y{x,t) and 
z{x, t) in a rotating Cartesian frame, where x G [0, /] is the spatial coordinate, 
I is the length of the beam. The above Cartesian frame is obtained from the 
fixed one by subsequent rotations on the angle (prit) (turning angle) and 
^Pnit) (raising angle). The system is controlled by torques Mj- and Mji 
applied at the bottom part of the beam. For each constant value ip^^, there 
is the control torque Mr = implementing an equilibrium (prit) = 0, 
^R(t) = ip% y{x,t) = 0, and z{x,t) = zo{x). The linearized system of 
differential equation describing oscillations around the equilibrium can be 
written as follows (see 



(1) 



1 / \ // 

y{x, t) + - [c;,y"{x, t)j = iIjt{,x)(pt, x G (0, /), 



1 / \ " 
z{x,t) + -\Cyz"{x,t)j = g(pRsmLp%^- 




(2) 
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—{czy")' - y + i^T{x) 
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= 0, 
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-c^y" - J^y' + Jj.i^'rpix 




= 
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— {cyz")' + QifR sin (p%- z- 




x=l 


= 0, 


(6) 


CyZ" + J2{^R + Z) 


x=l 


= 0, 




(7) 



{/o + (/i + Ji)sin2 ip^R + ih + h)cos^ mo{R-dcos^l.f+m{R-l cos^°Rf+ 
+ J {R — X cos ip^j^Y p dx^ ipT + J {R — X cos ip^j^)ypdx+ 
+ { mRy -{mly + .^y) cos 



Mt, 



(8) 



I2 + J2 + ^od + ml + / X pdx > ipR+ / zxp dx + < mlz + J2Z 
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zp dx + mz\^^i + ^^od + ml + J xpdx^ (pR^sirKf^j^— 

~ 9 zopdx + mzo{l)^ if R cos (p°ji = Mr, (9) 

where z{x,t) = z{x,t) - zo{x), (^R{t) = (^R^t) - Lp%, Mr = Mr- Mj^, and 

iPt{x) = X COS if^j^ — zo{x) sin (/j^j — R- 

We use dots to denote derivatives with respect to time t, and primes to 
denote derivatives with respect to the space variable x. The procedure for 
computing Zo{x) and M^^ is given in [T3] . 

The parameters in ([I])-© have the following physical meaning: p{x) is 
the mass per unit length of the beam, Cz{x) = E{x)Iz{x), Cy{x) = E{x)Iy{x), 
E{x) is Young's modulus, Iz{x) and Iy{x) are moments of inertia of the cross 
section of the beam with respect to the axes z and y, m is the payload mass, 
Ji, J2, and J3 are central moments of inertia of the payload, R is the platform 
radius, Iq is the moment of inertia of the platform, Ji, I2, and J3 are moments 
of inertia of the hub, mo is the hub mass, d is the distance between the origin 
of the rotating Cartesian frame and the hub center of mass. 

To simplify these equations we substitute expressions ([I]), ([2]), ©-([Tj) for 

y{x,t), z{x,t), y,y',z,z into ([8]), ([9]) and perform integration by parts 

x=l 

with regard for the boundary conditions ([3]). As a result, equations ([8]) and 
iQ take the following form: 

(fT = Ut, (pR = Ur, (10) 

where 

Ut = {Iq + {h + Ji) sin^ + mo{R - dcos(p'j^y+ 
+ (/3 cosv^^j + Jsz'oil) sinv3?j) cosv??^ + (^m{l cos ip^j^ - R)zo{l) + 

+j {xcosip°R-R)zopdx^ sin ip°R}-^x{MT- (i?(c,?/")' + c./ cos (/??j) , 

(11) 

Ur = {h+mod'^}^^ X ^^MR+Cy z'l^^^+g (^ j zpdx + m z\^^i + mod^ suiLp^^+ 
+ 9 (^j zopdx + mzo{l)^ LpRCosip^^f^y (12) 
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For each (puit), y{-,t), z{-,t), formulae (fTT]) and (fT^ establish a one-to- 
one correspondence between the torques (M^, Mji) and angular accelerations 
{ut, ur). Thus, we may consider [ut-, ur) G as a new control for the linear 
system (H])-©, (UnD- 



3 Main Results 

Consider the following linear space 

/ ^(0 \ 
C(-) 



X 



(f),uj,p,q GM. 



P 

V Q ) 

where H^{0, 1) is the Sobolev space of all functions whose generalized deriva- 
tives of order j = 0,1, ...,k exist and belong to L2{0, 1). For 

/ ??1 \ / ^2 \ 

'^^ ex and 6 = 

Pi P2 
\ Ql J \ <}2 / 

the inner product in X is defined by the formula 
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ex. 



X 



{ri'l{x)r]l{x) + C,i{x)C,2{x)) dx + 0i02 + UJ1UJ2 + P1P2 + qiq2- 



It is easy to check that the norm = ^/WTOx equivalent to the 

standard norm in H'^{0,1) x ^2(0,/) x (see, e.g.,P Ch. 3]), and hence, 
(X, II ■ \\x) is a Hilbert space. 

In order to consider an abstract formulation of the boundary value prob- 
lem ([I])-©, (ITOl) . let us introduce the linear operator A : D{A) X and the 
element B e X as follows: 

/ ^ \ / C \ / \ 



c 

to 
P 

V 1 / 



c 

i(cr^")" + 70 




B 





1 



(13) 
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where the domain of definition of A is 



C(o) = c'(0) = 0, }, 



(14) 



functions c(x) > and ip{x) are assumed to be of class C^[0, /]; J > and 7 
are constants. 

Let (^y{x, t), z{x, t), ifrit), j be a classical solution of the boundary- 
value problem ([I])-([7]), ([lU]) with controls {uT{t),UR{t)) for < t < r, r < 
+00. Defining 



( y{;t) \ 
iprit) 

ipTit) 

\ y'ii,t) I 



( 



\ 



~z{;t) 

'zil,t) 
z{l,t) j 



(15) 



we see that ix^t) e D{A) and ^/j(t) e B{A) for all t G [0,r). Consider the 
pair {At^Bt) obtained from (A, -B) by placing '^{x) = i'rix), c{x) = c^(x), 
= 7 = in f|T3|) . Similarly, let the pair {Ar,Br) be obtained from 
{A,B) by plugging ipi^x) = —x, c{x) = Cy{x), J = J2, 1 = gsimp^. Then 
the boundary- value problem ([I])-([7]), ffTOl) is reduced to the following control 
system: 

= At^t + BtUt, (16) 

in = Ar^r + BnUR, (17) 

where (^t, ^r) is the state and {ut, ur) is the control. In the sequel, we treat 
this control system as an abstract formulation of ([I])-([7]), (fTOj) with ^r, ^i? G X 
and Mr, ur G M. We see that (ITBl) . (IT7I) is separated into two parts, therefore, 
the stabilization problem may be solved independently for and ^r. The 
basic result we shall prove is the following 

Theorem 1. Consider the abstract Cauchy problem on t > 0: 



i{t) = Am + Bu, 



(18) 
(19) 
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where A, B are given by Ip^) . 



u = -— {ku + { a - 'J 



pip dx + mip{l) J 1 0+ 



+ f cr]"ilj"dx + [cr]"ip' - {cr]")'iP^ 



z=0 



- 7 



prjdx + mri{l) 



(20) 



a > 0, f3 > are large enough constants, and k > is arbitrary. 

Then the Cauchy problem [W\) . ( fTPj) with l[2U\) is well-posed on t > (in 
the sense of mild solutions), and the feedback control ^MW strongly stabilizes 
the equilibrium ^ = of the control system l[T8\) . i.e., for every e > 0, there 
exists 6 = 6{e) > such that, for every solution of (I^-^2^), 

\Mx<S^\\m\\x<e,\/t>0. 

Moreover, if a semitrajectory {^(t)}t>o of ( fl^) . is precompact in X then 
the set of its limit points (as t ^ +oo) is an invariant subset of Zq = {(^ G 
X : w = 0}. 

Proof. Consider a quadratic functional on X: 

2V{i) = + piu^+ [ I (C - i^tufp + v"^c] dx+ 



+ m{p-ip{l)uj}^ + J{q-ilj'{l)uj}^ -2-f(f)I^J rjpdx + mr]{l)j . (21) 
Let us compute the time-derivative of V along trajectories of (fT8|) when 

V{0 = {VV{0,A^ + Bu) = 

-I 

+ 



cC"v" - C ■ {cv")")dx + (p {cr]"y - qcv"^ 

+ + (3u + [cT]")" - p7(0V^ + r/) j dx+ 

+ {^ip'r]" — ip ■ {cT]")' — ni'-){(pip + 7])^ 



(22) 

By performing integration by parts with regard for conditions f lT^ . we get: 

Jo 
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c ■ {cv'y - Ccv") 



x=l 

I pI 

ijj ■ {cri")"dx = tp ■ {cr]"y\^^Q — / 'ip'{cri")'dx 
Jo 

I 



("cri"dx, 



■ {CT]")' - i^'cT]" 



x=0 



il)"cri'^dx. 







Let us substitute these formulae into fl22]) and use boundary conditions p = 
C{1), q = C'(0 from f[T^ . As a result, the expression for V takes the following 
form: ^ 



x=Q 



7 I / fyq dx + mrjil) ) >uj. (23) 



+ j cip"ri"dx + i^ip'rj" — ijj ■ {crj")'^ 

If u is defined by ( l20l) then formula ( l23l) yields 

= -kuj^ < 0, (fc = const > 0). (24) 

The next step is to prove that V{^) satisfies the following estimates 

Milieus <2V(0<M2||e||^ (25) 

with some constants < Mi < M2 < +00. On one hand, by exploiting 
inequahties (a + 6)^ < 2a^ + 26^ and 2ab <a? + b'^ in (pT]) . we obtain 

21/(0 < + ^"^^ + / + + ^'u;')p)rfx+ 



+ 7^0^ + 2(^1^ TIP c/x^ + 2 {mri{l)f . 
Then the Cauchy-Schwartz inequality implies 



r^p I < / rf dx I dx, 
J Jo Jo 



(26) 



(27) 
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r^2(/) = ^j^ ^'^a;^ < j^dx 7]'^dx. (28) 

The functions r]{x) and ri'{x) subject to the boundary conditions ri{0) = 
ri'{0) = satisfy Friedrichs' inequahties of the following form (cf. pET, p. 440]): 

I j2 f-l ^4 f-l 

Tfdx <- ri'^dx < - ri"'^dx. (29) 



2 Jo " 4 Jo 

By using inequalities fl27j) - fl29l) we conclude that 



i \ 2 rl 



rjp dxj + m 1] (l) < J T] dx J p dx+ 

+ Im^ r^'^dx < ^ + ^ P'^^) V"^dx. (30) 
Application of this inequality in (!26|) yields an estimate 2V{^) < M2\\C,T\\xy 
M2 = max< a + 7^, 2m, 2 J, 2 max p(x), 

(3 + 2 f ij^pdx + 2Ji:'\l) +2mij\l), 
Jo 

\ rn? H — [ p^c/x I + max c(x)l. 

On the other hand, we see that the inequality a? = (a—b+b)^ < 2{a—bY + 2b^ 
implies (a — 6)^ > a? /2 — b"^. By using the latter together with —2ab > 
-x^a^ - 67x2 (x ^ 0) in ([21]), we get: 

2^(0 > + (3io' + J^[cr]"' + - P^V)(ia;+ 



2 2/2 
-X 7 

X' 



— ( / ?7p(ix + mr7(/) ) > 
> (a - x^^) 02 ^ !^p2 ^:Lq^^^j^ C^pdx+ 



+ /mine -(m^ + - [ p^dx\ \ [ 'n"^dx. (31) 

We have also used the inequahty (l30l) here. From (13T!) we conclude that 
2V^(0>Mi||el||and 



, ^ . r 2 2 m J 1 . 
Ml = mm< a — X 7 , — , — , - mm p[x 
I 2 2 2 xg[o,«i 



(3- / ptp^dx-mtp\l)-J^'\l), 





mm c[X) I m 



provided that 



> — ^ -— (m^ + t: [ P^dx 1 , 



min^.g[o,i] c(x) V 2 







a>xV, /3> / pip'dx + mil)\l) + Jip'^l). 







For the rest of this paper, we assume that constants a, /3, and x satisfy the 
above inequahties. 

The estimate (125!) shows that the two norms and ||^|| y = a/V^(0 are 
equivalent in X. Let us write the closed-loop system (fT8|) with the control 
M defined by (!20l) as ^ = A^, where -D(^) = -D(A) is dense in X. From 
inequality fl21|) it follows that the operator A is dissipative in X equipped with 
the norm || ■ ||y. Then the Lumer-Phillips theorem [21 Chap. 1.4] implies that 
A is the infinitesimal generator of a Co semigroup of contractions, {e*^}(>o, on 
X (with respect to the norm || ■ ||y). It means that the Cauchy problem (fT8|) - 
( 120|) has the unique mild solution ^(t) = e'^^^o, t > 0, for every ^ and 
the above solution is classical if ^ D{A). As {e*^}t>o is contractive (under 
an equivalent renormalization in X), then 

\m\\v<mv. vt>o. 
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This implies, taking into account the estimate that 



The above inequahty proves strong stabihty of the equihbrium ^ = in the 
sense of Lyapunov (we may choose 6{e) = e a/ Mi / M2 in the definition of 
stabihty). 

To conclude the proof we apply LaSalle's invariance principle |S1 [ID] 
(cf. [m Lemma 2]) with the functional V{^): if a semitrajectory {^(t)}t>o 
is precompact then its cu-limit set, fl{C,o), is a non-empty and semi-invariant 
subset of G D{A) : V{^) = 0} = Zq. □ 

Remark. By applying the formula fl2Ul) to control systems ([TB]) . flTTl) 
separately and using the representation (fT5|) . one can write the feedback 
control proposed as follows: 



x=0 



CyZ 



\x=0 



+g { / {xifn — z)pdx + ■m{l(pR — z\ 



x=l. 



+ 



sin 



To implement these controls in practice, it is sufficient to compute ut and ur 
depending on the measurements of ipx, ^Pr, (fr, <fR, y, z at each t > 0, and 
then apply formulae (fTTj) . (fT2!) to find torques Mt and Mj^. An advantage 
of this approach is that no information about the time-derivatives of y{x,t) 
and z{x,t) is needed. 



3 Conclusions 

A feedback control has been derived to stabilize the equilibrium of a dif- 
ferential equation in a Hilbert space that describes the motion of a fiexible 
beam with a tip mass. Although the main result of this paper concerns 
non-asymptotic stability, further analysis of the limit behavior of controlled 
trajectories is possible by means of the invariance principle. The main diffi- 
culty in this direction is to prove that the semitrajectories are precompact. 
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which is not an easy task in general (see, e.g., [21 We do not study the 
compactness issue here, leaving it for future work. 

This research is supported by the Ministry of Education and Science of 
Ukraine through grant of the President of Ukraine for Young Scientists. 
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